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Abstract. We consider d — 2 Ising strip with surface fields acting on boundary spins. 
Using the properties of the transfer matrix spectrum we identify two pseudotransition 
temperatures and show that they satisfy similar scaling relations as expected for 
real transition temperatures in strips with d > 2. The solvation force between the 
boundaries of the strip is analysed as a function of temperature, surface fields and the 
width of the strip. For large widths the solvation force can be described by scaling 
functions in three different regimes: in the vicinity of the critical wetting temperature 
of 2D semi-infinite system, in the vicinity of the bulk critical temperature, and in the 
regime of weak surface fields where the critical wetting temperature tends towards 
the bulk critical temperature. The properties of the relevant scaling functions are 
discussed. 



PACS numbers: 05.50.+q, 68.35.Rh, 68.08.Bc 
1. Introduction 

Fluctuating condensed-matter systems enclosed by walls are characterised by the 
appearance of solvation force acting between the walls. This force originates from the 
fluctuations of the confined system. The properties of solvation forces have been the 
subject of increasing interest during the last years [Tl[2l[3llll[5l[6l[71[8l[9l[l0l[ni[l2l 
[T3l [HI [151 [m [13 [IB]- Both the shapes and possible chemical inhomogeneity of the 
confining walls influence the form of solvation forces [H M, [IS, [161 IE] which additionally 
depend on the thermodynamic state of the system and on the interaction between the 
system and the walls. In particular, if the system is chosen to be at its bulk critical 
point the solvation forces become long ranged and show universality [19] while in the 
vicinity of criticality scaling behaviour is observed [20l [21] . 

In this article we analyse the solvation forces in two-dimensional Ising strips. The 
spins are confined by two parallel, planar and chemically homogeneous walls separated 
by distance M. Each wall interacts with the system by surface fields acting on the 
boundary spins. Our goal is the exact determination of the properties of the solvation 
forces as functions of temperature, surface fields and the width of the strip M. To 
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investigate these properties we a use method based on exact diagonahzation of the 
transfer matrix which is followed by numerical determination of appropriate eigenvalues. 

The paper is organised as follows. In section [2] we define the 2D Ising strip, recall 
its properties in the bulk limit as well as the properties of the semi-infinite system 
related to critical wetting. In section [3] we define pseudotransition temperatures and 
check that they display the scaling properties expected for higher dimensional systems. 
Section H] is devoted to our main goal, i.e. analysis of the properties of the solvation 
force acting between the system boundaries. We first recall the definition of solvation 
force and adapt it to our model. We study this force numerically to establish several of 
its properties as functions of temperature, surface fields and the distance between the 
walls. For large width of the strip we explain these properties by introducing scaling 
functions in three different scaling regimes: around the wetting temperature, around 
the bulk critical temperature, and in regime in which both of the above temperatures 
are close to each other. 

2. Ising strip 

2.1. The model 

We consider an Ising model on a two-dimensional square lattice with columns and 
M rows, and impose periodic boundary conditions in the horizontal direction. In this 
way the Ising strip of width M is obtained. We assume that surface fields hi and h2 
act on the spins located at the bottom and the top row, respectively; these fields can be 
considered as model short range interactions between the system and the surrounding 
walls. The Hamiltonian of the system has the form 

N M-1 N 

({Sn,m}) = ~J {Sn,mSn+l,m + Sn,mSn,m+l) ~ ^^(^I'Sn,! + h2Sn,M)y (l) 

n=l m=l n=l 

where Sn,m = ±1 denotes the spin located in the n-th column and m-th row, and 
SN+i,m = Si^m- The coupliug constant J is positive (ferromagnetic case) and we assume 
no bulk field h acting on the system. 

In this paper, we concentrate on two special choices of surface fields corresponding 
to the so-called symmetric and antisymmetric case: in the symmetric case (denoted by 
the superscript S) one has hi = h2 while in the antisymmetric case (AS) hi = —h2. 
Later on, we will also use superscript O to denote the limiting case hi = h2 = which 
is referred to as the free case. 

2.2. The free energy of the strip 

To calculate the free energy of our system we use the method based on exact 
diagonahzation of the transfer matrix. In this method the 2^^"*"^ x 2^^"*"^ transfer matrix 
is represented by (2M + 4) x (2M + 4) orthogonal matrix R. Eigenvalues of the transfer 
matrix are calculated from eigenvalues of R which can be found by solving recurrence 
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equations for eigenvectors of R [221 ESI 121] • Here we only recall the final formulae for the 
free energy per column (here and in the following formulae we do not explicitly write 
the dependence of the free energy and other quantities on the coupling constant J) 



f{T,h,M) 



-kuT 



1 M 
- (7i + 72 + 73 + • • • + 7m+i) + yln (2sinh2/s:) 

1 M 

2 (-71 + 72 + 73 + • • • + lA4+i) + Y 111 (2 sinh 2K) 



(2) 

(3) 
. < 



where K = J/k^T and k-Q is the Boltzmann constant. The coefficients 7i < 72 < 
7Af+i are positive functions of parameters T, hi and M defined by relation 

cosh7j!c = cosh {2K — 2K*) + 1 — cosuk, (4) 

where parameter K* is obtained from sinh sinh 2i^* = 1. The functions Uk are 
solutions of the equations 

{M + l)uk-6'{uJk,T)-(f){uJk,T,hi) = {k-l)7r, < cu^ < vr, (5) 

and k = 1,2, ...,M + 1. The function / (T, hi) is defined as 



l{T,hi 



2 forr<Tw, 

for < T < Tr. 



1 for T>Tc. 

The symbol denotes the bulk critical temperature [25] 

J __ 
~ ~ 2 



(6) 



(7) 



while Tw {hi) denotes the temperature of the critical wetting transition taking place in 
the semi-infinite Ising model. It depends on the surface field hi and can be defined by 
equation [26] 

W{T^,hi) = l, (8) 

where 

hi 



W (T, hi) = (cosh 2K* + 1) (cosh 2K - cosh 2Ki) 



knT' 



(9) 



We observe that for certain ranges of temperatures ([5]) may not have a solution for k = 1 
and k = 2. In such cases ui and 002 are imaginary and satisfy equations 



LOk = lUk, 



-UkM 



ak exp {i [(f) (iuk, T) + S' {iuk, T, hi)]}, 



1,2, 



(10) 



with ak = ±1. The detailed rules for selecting the signs of ai and 0^2 are presented in 
the next subsection. Functions 6 and 6' are calculated from the formulae 



, - 1 



(e''^ - A) (e''^ - B) 



:iii 



(Ae^^ - 1) {Be^ - 1) ' 

where A{T) = (tanhiT tanhi^'*)"\ B {T) = tanhiT/tanhi^'*, and the function 
W{T,hi) is given in ([9]). To determine the angles (f){uj,T,hi) and 6' {uj,T) uniquely 
we pick the continuous branches of solutions for which 

0(O,T,/ii) = 7r, 5'(0,T) = -7r. (12) 
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For u = and T = T^, the angle (j){uj,T,hi) is undefined (W (T„, hi) = 1) while 
at T = Tc the angle 6' {uj,T) is undefined {K = K*, so B = 1). Although at these 
temperatures our formulae are useless one can use the continuity of the free energy and 
calculate it using a limiting procedure. 



2.3. The characteristic temperatures 

Because the lower critical dimension of the Ising model equals two [di = 2), no true 
transition may occur in a two-dimensional Ising strip with finite M. On the other hand, 
the infinite 2D Ising model experiences the critical point behaviour at T = Tc, while in 
the semi- infinite 2D Ising model with the surface field hi, the critical wetting transition 
takes place at T = T^, < Tc [20\. Below, we discuss the properties of ([5]) and on this 
basis we define the characteristic temperatures Tjj^j and Tj^. 

First we consider T < T„ case, for which 1 = 2, see ([6]). For small enough 
temperatures the left-hand side of ([5]) is an increasing function of u, it equals for 
uj = 0, and thus this equation does not have a solution for k = 1 and k = 2. The 
coefficients cui and 002 are thus found from (|TOl) with ai = —1 and 02 = 1- However, 
when T is getting close to T^ the situation becomes different: the left-hand side of ([5]) 
— upon increasing u — first decreases, has a minimum and then increases. As a result 
([5]) has a solution for = 2. At the same time, to obtain coefficient ui equation ( ITOj) 
must be used with ai = —1. The M-dependent temperature, which separates the above 
two possibilities is denoted by Tj^. 

When T^ < T < Tc, one has / = and all coefficients Uk are defined by ([5]). 

For T > Tc, I = 1 and for temperatures well above Tc equation ([5]) does not have 
solution for = 1; the coefficient Ui can be calculated from (fTOj) with ai = +1. When 
T is close to Tc the left-hand side of ([5]) is a non-monotonic function of u, and thus the 
solution exists for any k. The characteristic temperature separating these two cases is 
denoted by T^j^,j- 

Typical plots of the left-hand side of ([5]) are shown in figure [H 

To find the formulae for M-dependent temperatures and Tj^^ we use the fact 
that at these two temperatures ([5]) has double solution for u = 0. In other words, the 
condition 

d_ 

duo 

must be satisfied, which leads to 

2W{T,hi) sinh2ir ^ 

^ ' = M + 1, f 14) 

W{T,hi)-l smh{2K-2K*) ' ^ ' 

where W (T, hi) is defined in ([9]). To find solutions of this equation it is useful to analyse 
its left-hand side as a function of temperature: it equals to 1 for T = 0, is an increasing 
function of temperature for < T < T„, at T = T„ reaches infinity and has a pole 
{W = 1 for T = T„). For T^ < T < Tc the left-hand side of f|T^ is negative and has 
another pole for T = Tc {K = K* at T = Tc). For T > Tc it decreases from infinity at 



[(M +1)lj-6' {u, T)-(f) {u, T, hi)] = (13) 

uj=0 
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Figure 1. Plots of the left-hand side of (O divided by tt for hi = 0.8 J (Tw « 0.621 Tc) 
and M — 5. The solutions ujk correspond to integer values of this function. For 
T = O.STc < T^^j^i there are A/ - 1 = 4 solutions, for Tj_^^ < T = 0.6Tc < T„ there 
are M = 5 solutions, for < T = 0.8Tc < Tc there are M + 1 = 6 solutions, for 
Tc<T= l.OlTc < T^'j^ there are M+1 = 6 solutions, and for <T= 1.2Tc there 
are M = 5 solutions. Total number of solutions is A/ + 1 = 6; the missing solutions 
correspond to imaginary values of lu and are determined from (|10p . 

100 I , n. , , , n , 1 
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Figure 2. Plot of the left-hand side of ^ for hi ^ 0.8 J (for which Tw « 0.621 Tc). At 
pseudotransition temperatures {hi,M) and j^j {hi,M) this function is equal to 
M + 1 . To guide an eye the vertical broken lines corresponding to T = and T = Tc 
are drawn. 

T = Tc to for T ^ oo. A typical plot of left-hand side of f[T^ is shown in figure [2J 
Equation (fT4|) has two solutions for any positive M - the solution Tj^^ is always smaller 
that T„ and approaches the wetting temperature monotonically as M — >■ cx), while the 
solution T^j^ifii^M) is always larger than and decreases monotonically to Tc as 
M ^ oo. 

3. Properties of pseudotransition temperatures 

In an infinite strip of width M and dimension d larger than the lower critical dimension 
di, d > di = 2, true phase transitions corresponding to the non-analyticity of free energy 
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Figure 3. The plot of the amphtude Ai (see (fT5|) ) as a function of 



occur. 

In a strip with symmetric surface fields (S) capillary condensation is expected. For 
a vanishing bulk field the strip is filled with phase favoured by the walls for T < Tc^m- 
The critical temperature T^^m is shifted away from T^. |27j. On the other hand, for 
antisymmetric surface fields (AS) with no bulk field, a transition is observed at T^^m 
that is shifted from T^. For T < Tw,m the interface separating two phases is located 
close to one of the walls while for T > T^^m this interface is located in the middle of 
the system [20]. Temperature T^^m approaches as M — > oo. In AS case the second 
phase transition at Tc^m located close to Tc, also occurs. 

In a two-dimensional strip no phase transition may occur for finite M. However, 
for large strip widths we expect some thermodynamics functions to vary rapidly 
close to certain temperature values while remaining analytic. It is convenient to 
define these pseudotransition temperatures which can be then used to characterise the 
behaviour of our system. Since all functions are analytic, these temperatures cannot be 
defined uniquely. There are different criteria according to which the pseudotransition 
temperature can be defined and thus there is no single T^^m and T^^m- One possibility 
corresponds to and T^m defined as the temperatures at which the specific 

heat attains its maximum values. Here we would like to show that the just defined 
temperatures and Tj^j may be treated as such pseudotransition temperatures. 

First we check how the difference T„ — Tj depends on the width of the strip M 
for large M. This can be done on the basis of f|T4l) . Using the implicit function theorem 
one obtains 

(15) 



M 



T=T„ 

Figure [3] shows the plot of the amplitude Ai (T„) after reparametrization from hi to T„ 
has been done according to ([8]). 

Parry and Evans [20j used scaling hypothesis to postulate that ioi M ^ oo 
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Because for a 2D Ising model (3s = 1 the behaviour of the difference between 
pseudotransition temperature and agrees with this hypothesis. 
Similarly, for one obtains from (HM 

= ^ + O (M-) , (17) 

where the amplitude A2 = [2 In (l + V^)] ^ is universal. Since for a 2D Ising model one 
has 1^ = 1, thus {T^m — Tc) /Tc ~ M'^^'^ as expected on the basis of scaling arguments 
[2T] . We note that Tj^^ is always larger than Tc. 

The wetting temperature is a continuous function of the surface field hi. Parry and 
Evans [20j proposed the scaling function Xas which describes the dependence of Tw,m 
on the width of the strip M and the surface field hi in the limit /ii — >■ and M — 00 
with /iiM^i/'^ fixed 

~ ^"''^ = M-^/'-Xas {hiM"^'/^) . (18) 

^ c 

It turns out that the pseudocritical temperature ^ defined by ffT^ satisfies a similar 
scaling relation. We have found the exact expression for the corresponding scaling 
function X^g. 

For a 2D Ising model Ai = | and the scaled variable takes the form x = /iiM^/^. 
In order to find the scaling function XTg (x) 



^ = M-'Xl, (x) + O (M-^) , (19) 

c 

we introduced in (IT^ the surface field hi = xM"^/^ and obtained in the scaling limit 
X2g(x)= [21n(l + v^)]"' + i(l + V2)ln(l + v^) (20) 

It is interesting to note that this result is true only for fixed x; the scaling function 

Tc - TJ (M, xM~i/2) 
lim M- "'^^^ ' ^ (21) 



^7 / 



M~^oo T 



c 



is not a uniform limit. 

The scaling law ( JT9l) has finite size corrections of order which are present even 
for hi = 0. 



4. Solvation forces 
4.I. Definition 

The free energy of the strip per column can be calculated from ([2]) and ([3]). For both S 
and AS cases it naturally decomposes into the sum of three terms 

r (T, hi,M) = Mfi, (T) + f: (T, hi) + /f,, (T, /ii, M) , (22) 
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where a G {S, AS}, /b is the bulk free energy density [2S! equal to 

(T) = -ksT — / arccosh [cosh {2K - 2K*) + 1 - cosu;] do; 
L27r Jo 

+^ln(2sinh2X) , (23) 

/g" (T, hi) is the surface free energy per column, and the remaining term f^^ (T, hi, M) 
describes the interaction between the boundaries of the strip per column. By definition, 
the surface free energy f^ (T, hi) does not depend on M, and f^^ (T, hi, M) tends to 
as M — i> oo. 

In general, the solvation force is defined as minus derivative of /"^ with respect to 
the distance between the boundary walls. In the present case, because M is integer, we 
use the definition 

(T, hi, M) = - [/;-,, (T, hi,M + l)- fr^, (T, hi, M)] /k^T, (24) 

where the factor l/k-^T is additionally introduced to make the solvation force 
dimensionless. This definition is equivalent to 

{T, hi, M) = [r (T, hi,M)- r {T, hi,M+l) + f, (T)] /ki,T. (25) 

It is also useful to introduce the difference between the solvation forces 
corresponding to different boundary fields configurations 

A/,oiv (T, hi, M) = f^^ (T, hi, M) - /^i, (T, hi,M). (26) 
Using ([2]), ([3]) and ( l25l) it is straightforward to show that 

A/sol. (T, hi, M) = 71 (T, hi,M)- 71 (T, hi,M+l). (27) 
This difference is easier to study analytically than the expression for /^^j^ (T, hi, M), see 

m. 



4-2. Basic properties 

We start our analysis by evaluating numerically the solvation forces for different 
temperatures T, strip widths M and surface fields hi. 

In the symmetric case (S) the solvation force is always negative (attractive). For 
hi close to J this force has a minimum at T^f„ > Tr and tends to both in the small 
and large temperature limits. Upon decreasing the boundary field hi, the absolute 
value of solvation force decreases, and for hi small enough a second minimum appears 
at T^J^ < Tc. Upon further decreasing hi, the minimum located at disappears. 
The range of hi for which /^^j^ has two minima depends on M, and for M ^ oo this 
range shrinks to 0. Plots of the solvation force in the symmetric C3iSG clS cl function of 
temperature for different boundary fields are presented in figure |H The behaviour of 
this force will be studied in detail using scaling functions later on. 

In the antisymmetric case (AS) the solvation force is plotted in figure [51 For hi = J 
this force is positive (repulsive) for all temperatures and has maximum at T^^^ located 





Figure 5. Plots of the solvation force in the antisymmetric case (7i2 = —hi) as a 
function of temperature for M — 25 and different values of the boundary field hi . 

slightly below Tc. The solvation force f^f^ (T, hi = J, M) tends to in the high and 
low temperature limits. However, for hi < J the solvation force changes sign. It is 
negative for small temperatures, has zero at T* slightly 

above T^. For temperatures higher than T* the solvation force is positive and has a 
maximum close to T^. For hi approaching 0, T* tends to and the (negative) value 
at the minimum below the wetting temperature decreases. The (positive) maximum 
value of the solvation force also decreases and disappears in the limit hi 0. We also 
looked at the location of the maximum of the solvation force T^£,. For small M, T^^^ 
is located above T^. Upon increasing M the temperature first crosses the critical 
temperature and then, upon further increasing of M, approaches T(. from below. The 
exact value of M at which Tj^f, is equal to Tc depends on the boundary field hi. We note 
that the limiting value of the solvation force at /ii = is the same for both boundary 
fields configurations. Some of the above described properties of solvation force have 
been reported for a different system in [9]. 
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The leading M-dependence of the solvation force evaluated at is known exactly 

IT 



/i, {T,,h,M) 
/o, (T„0,M) 



48M2 
48M2 ^ W J 



+ O (l/M^) . 



(28 a) 
(285) 
(28 c) 



(30) 



48M2 

The above values of the universal amplitudes are also recovered numerically in our 
analysis. 

We checked numerically that for T ^ [31] 

fl,^ (T, h, M) ~ exp [-M/^b (T)] , (29) 

where [32] 

UAK-AK*y^ forr<r, 
\ (2ir* - 2Ky^ for T > T, 

is the bulk correlation length. Using fl271) and the dependence of 71 on M for fixed T 
[26] we checked that ( 129|) implies the following leading order decay of the solvation force 
in the antisymmetric case 

'exp [-MlnPF (T, hi)] 
exp [-M/i^{T)] 

^exp [-M/eb(r)] 

The solvation force is a continuous function of temperature and the above formula is 
correct only in the M — > cxd limit. Below we discuss the behaviour of the solvation force 
around T„ and Tc by introducing the appropriate scaling functions. 



f^^{TM,M) 



for T < Tw, 
for T = T^, 
for T„<T <T^ 
for T = Tc, 
for T > Tr. 



(31) 



^.5. Scaling at T„ 

To study properties of the solvation force close to T„ in the antisymmetric case we take 
the scaling limit M —>■ 00, T ^ with parameter X = MlnW (T, hi) ~ (T„ — T) M 
fixed. The function W (T, /ii) has been introduced in the scaling variable to simplify 
the scaling function. 

To study the solvation force in this limit we use (127|) . For T < T„ coefficient 71 is 
given by (fTOi) with k = 1 and ai = — 1 

(T, /ii, M) e-" (e-" - W (T, /ii, M)-^) 



-uM 



where the solution u gives 71 using 



e-u -W{T, hi,M) 
with uJi = iu. 



(32) 
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We put M = X/ \nW in (15^ and calculate the limit T ^ using I'Hopital's rule. 
After introducing H (X) = {W, X)\^__^ one gets 

e-^^W = |||l^, HiXXO. (33) 

The function H {X) can be calculated numerically for any X > 0. The solution of (!32l) 
in the scaling limit takes the following form 

u = H {X)\nW + {in^W) , (34) 

and one obtain^ 

(X, M) = . (/,) - ^^1^ + O (^) . (35) 

where !/ (fti) = (2A--2A-)|r,j,^,j,,. 

Comparison of this result with fl29|) leads to the conclusion that the first term on 
the right-hand side of ( l26l) dominates in the scaling regime (except of T = T^) and the 
second term may be neglected. 

With the help of (135|) one gets 



where 



/stfv (X, h,M) = -—G (X, hi) + Ol^—), (36) 
C(Xh)- ^"^'^^^ g^(X)-l 

^^^'^^^-sinhz.(M 2 + XiHHX)-iy ^^^^ 

Using (!33l) and (!37|) it is straightforward to analyse the properties of the scaling 
function G (X, hi). For small X 

G(X,/ii) = -^ + 0(x2), (38) 
smh z/ ^ 

it has a maximum at Xq ~ 3.22149 and approaches zero exponentially for large X, see 
figure [61 

The above properties can be used to explain the behaviour of the solvation force 
around the wetting temperature for large M. In particular, one has 



fiS, (r- . k, M) = -MM + o (-Ij) , (40) 

where the functions ^3 {hi) and ^4 {hi) are positive and may be obtained from the 
scaling function G and the definition of scaling variable X. 

The dependence of T* on M can be explained using (!26|) . Exactly at T = the 
left-hand side of this equation is zero (exactly at T„ the coefficient 71 {T,hi,M) does 

I Note that the expression p5)) for the function 71 (X, /ii,M) is not equivalent to equation (4.7) in 
[24] . because of an error in calculation. 



Properties of the solvation force of a two-dimensional Ising strip in scaling regimes 12 



2.5 I > < < < < < < < r 




0123456789 10 
X 

Figure 6. The scaling function G{X, hi) for the solvation force in the antisymmetric 
case multiplied by sinhi/(/ii) (see ([55)) ). 



not depend on M) so the solvation force is the same for both AS and S boundary fields. 
With the help of equations (ETj), (I29l) and (EH]) one gets 

^^^^ = A, ih) M' exp [-M/a (T„ (h))] , (41) 

c 

where {hi) is a positive function. 
4-4- Scaling at 

For temperatures close to the bulk critical temperature the solvation force takes in the 
limit T ^ Tc, M — > oo with fixed hi and x = sign (T — Tc) M/^;, (T), the following 
scaling form 

/.V (T, hi,M) = -^X^^ (x) + O (M-^) . (42) 

Note that the factor sign (T — T^.) introduced in the definition of the scaling variable x 
makes it negative for T < Tc and positive for T > T^. The bulk correlation length ( l30ll 
close to Tc takes the form 

|tr^ for T >Tc{x> 0), 
^oltf^ for T < Tc (x < 0), 

where t = {T — Tc) /Tc and the amplitudes = 2^^ = 1/ [21n (l + V^)] , such that 
X ~ (T — Tc) M for T close to Tc. Later on we will use 

fx for X > 0, 

X = tMlit ^{ ,^ , - n 
I x/2 for X < U, 

instead of the scaling variable x. 

The scahng function X has already been proposed by Evans and Stecki [31] and has 
been calculated analytically in both S and AS cases for particular value of the scaling 
field hi = J. Here we consider arbitrary values of hi. Our numerical calculations show 
that, up to numerical errors, 

X;:^ix)=X^ix) for hiy^O. (45) 



e.m^r._!' 1 . ^ J: S (43) 
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We have found numerically that the corrections to scaling depend on hi and are getting 
smaller for hi close to J. Unfortunately, we are unable to prove analytically this property 
of scaling function (x). 

However, the difference between the scaling functions for the two cases 

AX,, (x) = A-f (x) - (x) (46) 

can be calculated exactly. For T close to Tc this function is obtained from (127|) . (jlj), ([5]) 
and (E]). To derive uJi in the scaling limit we replace M with x^q /t in ([5]) and use the 
following property of the function (j) {u, T, hi) 

Urn (f){uJi{T,hi,M) ,T,hi) = (47) 

for hi 7^ 0. Thus the only term that depends on the surface field in ([5]) disappears in 
this scaling limit and the calculation of AX goes along the same lines as in [31] (from 
now on we drop the index hi in AX). This is in full agreement with (H5|) . One obtains 

^ ^, sinw , , 

AX = , 48 

w — sm w cos w 

with w being a solution of 

w cot w = X, (49) 

where 0<w<7rforx<l, and w = iu, u > for x > 1. 

Function X^_^ (x) for different /ii-values is plotted in figure [T^i. Note that such 
obtained curves are indistinguishable from each other which numerically proves for 
the symmetric case. It has a minimum for x > 0, so from fH2|) it follows that for large 
M 



Til {hi, M) =T, 



1 + ^ + (M-) 



(50) 



/soiv (r^fn, hi,M) = + , (51) 

with Aq and A-/ determined by the position of minimum of the scaling function 

Ae ^ 1.26424, ^ 0.43052. (52) 

Because X^^ (x) has only one minimum, the second minimum of the solvation force, 
located below Tc, disappears in this limit. 

Function X^^ (x) for for different /ii-values is plotted in figure [7b. Again one notes 
that such obtained curves are indistinguishable from each other which numerically proves 
( 145|) for the antisymmetric case. It has a maximum for x < and from fjl2|) it follows 
that for large M 



T^lihi,M) =Tc 



1-^ + (M-) 



(53) 



/stfv&/^i,M) =i| + 0(M-3), (54) 



A, 

with and Ag determined by the position of maximum of the scaling function 

^ 0.2651, Ag ^ 1.5341. (55) 
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Figure 7. The scaling function X describing the solvation force for any /ii ^ and 
Af — > oo with X — tM/^Q fixed for: (a) symmetric {hi — /12), and (b) antisymmetric 
{hi = — /12) boundary fields. Each plot does not depend on the chosen value of hi 
and is the same (up to numerical errors smaller than the resolution of the plot) as the 
analytically calculated scaling functions for hi = J [3T] . 



The temperature T^^^ is smaller than Tc in this limit. These results have already been 
reported in [31] for hi = J. According to our numerical analysis the values of constants 
Aq, Aj, Ag and Ag are the same for any nonzero surface field hi. [§| 



4-5. Scaling for T„ 



c 



To explain the properties of the solvation force for small values of the boundary field 
hi we consider the scaling limit M — 00, T — > Tc and /ii — (i.e. T„ — > Tc) with two 
scaling variables 

tM Ao hi 

fixed. In this limit the solvation force can be described by scaling function 3^" {x,y) 

(T, hi, M) = -^y- {x, y) + {M-') . (57) 

1 /2 

The constant = [(l + v^) / In (l + ^2)] ' in was chosen such that for negative 
values of x, the value y = I corresponds to T = T^. For y < I equation ( 1571) gives the 
solvation force for T below T^, and for y > 1 - for T above TL. This scaling function 
has already been analysed for subcritical temperatures in [18]1||1 

The scaling function y"' {x, y) can only be calculated numerically; details of 



evaluation are presented in Appendix A 



§ Note a minor disagreement between values of our numerical amplitudes Ag and Aq and those 
evaluated in [31j due to minor numerical inaccuracies in |31j . 

II There is a mistake in the scale of variable x in figures 3, 4 and 5 in [T^. To get the correct values of 
X one should replace x by [^q) x in these figures in [18]. 
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Table 1. The rules for picking the correct solutions of ([60|) corresponding to 71. In 
all four cases there exists exactly one solution for a given range of w. 



Range of x and y 




Domain of w 


Condition on w 








y > l,x < {y^ - 1) 


1 {v" + 1) 


real 


< w < niin ^tt, |x| | 






} 


y>l,x> {y^~l) 


/ if + 1) 


imaginary 


< w/i <x{f + 1) 


(1 + 


2yr 


-1/2 


y <l,x > 


imaginary 


X (1-2/4)1/2 < w/i 














w/i < a; (2/^ + 1) (1 ^ 








y <l,x<0 




imaginary 


|a;| (1 - 2/2) < w/i < 


\x\{\ 





Before presenting the numerically evaluated properties of the scaling functions 
y"(a;, y), a = S,AS we note that one can test some of these properties through 
analytically determined difference 

Ayix,y) = y^^x,y)-y'{x,y). (58) 

This can be done with the help of fl27|) . The coefficient 71 is given by equation ([5]), 
where its solution ui determines 71 by (jll). After applying the scaling limit to the above 
equation one gets 

71 (T, h,,M) = j^V^^T^ + O (M-2) , (59) 
where w is a solution of 

X2 (?/4 _ 1) _ ^2 y ^ 

Depending on x and y equation (!60l) may have many different solutions for w. The 
rules for choosing the correct solution are summarised in table [U other solutions give 
the coefficients 7^ for k > 1. 

The function A3^ is given by the formula 

which leads to a rather lengthy expression and we refrain from presenting it here. 

The scaling functions {x, y) are plotted in figure [HI These plots cannot be used 
directly to approximate the behaviour of the solvation force as a function of temperature 
for large fixed M, because for fixed y both the temperature T and the surface field hi 
become functions of x. Additionally, the limit a; — > corresponds to hi — > 0, which 
explains why — for any y — one has 

3^^'(0,y) = 3^^(0,y) = -^, (62) 

i.e. in this limit the scaling function equals the universal amplitude describing the decay 
of the solvation force at T = Tc for free boundary conditions. 

To explain the observed properties of solvation force we changed variables in the 
scaling function and defined new function 

3^"(x,z) = 3^"(a;,v/i7R), (63) 
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Figure 8. The scaling function y" {x,y) as function of x for different values of y. 
Boundary fields are symmetric (a = S) in graphs (a) and (b), and antisymmetric 
(a = AS) in graphs (c) and (d). 

where the new variable z = |x| ~ Mhl, so that fixing x and z is equivalent to fixing 
X and y in the scaling limit. For the new scaling variables one obtains 

(T, h, M) = -^y^ (x, + O (1/M3) , (64) 

which can be used to approximate the solvation force as a function of temperature for 
fixed M and hi. Plots of y are presented in figure [9l 
We checked that up to our numerical precision 

lim (x, z) = ;f ° (x) , lim3^"(x,z) = A'°(x), (65) 

>oo 2— >0 

which is not surprising since for M — > oo and fixed hi ^ Q one has z ~ Mh\ oo, 
while hi = implies z = 0. 

In both S and AS cases the scaling functions y°' {x, y) refiect the behaviour of 
solvation force for small surface fields. For symmetric surface fields, the negative function 
y^ {x, z) has for fixed z < zi 0.1474 only one minimum at negative values of x. At 
z = Zi the second minimum located at positive x appears. Upon further increasing of 
2;, the minimum at x < is increasing and the absolute value of the second minimum 
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Figure 9. The scaling function y" {x, y) for fixed z as a function of x. Boundary 
fields are symmetric (a = S) in graph (a) and antisymmetric [a = AS) in graph (b). 



is increasing. Wc could not observe the disappearance of the minimum at negative x, 
because for large z the depth of this minimum is of order of our numerical errors. 

For antisymmetric surface fields the scaling function (a;, z) has exactly one 
minimum and one maximum for all finite z. The minimum is always located at a; < 
and moves towards — oo when z is increased. The maximum moves from x — ooior z — Q 
to a finite negative value of x for 2; = oo. The value of scaling function at maximum 
is always positive and becomes very small for z close to 0. For z — Z2 ~ 0.212 the 
scaling function vanishes at a; = 0, which means that (up to higher order corrections) 
the solvation force disappears at T = Tc- On the other hand, at 2; = ^3 ~ 3.35, the 
maximum of scaling function is located exactly at x = 0. 

The above observations are summed up below: 



in the S case the minimum of the solvation force located above exists for 



xfM ^ ' 



0.40, 



in the AS case the solvation force is zero at T = Tp for 



Ill ~ 



0.48, 



in the AS case the maximum of the solvation force is located exactly at T 

Av2 



hi/ J 



+ 0(M-3/2), ^ 



12 



1.89. 



Finally, we study the behaviour of the solvation force at T 

{z) 



M2 



+ O (M-^) . 



The amplitude A"' {z) can be calculated using the scahng function 

[z) = y'^{x = 0, z) . 



(66) 

(67) 

Tc for 
(68) 

(69) 
(70) 
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Figure 10. The amplitudes A"" {z) describing the decay of the solvation force at 
T — Tc for M — > oo (see (|69| ): (a) symmetric surface fields, (b) antisymmetric surface 
fields. Dashed lines show the exactly known values of amplitudes for z = and z ^ oo. 



From (128 gp . (128 6|1 and (128 cp the values of the amphtudes ^4° (z) for z 
follow 



AS 



0) = A" {z 



OO = — - 



A 



z = oo] 



and z ^ oo 



(71) 



48 ' '48 

For other values of z amplitudes A°' (z) can be calculated numerically; they are presented 
in figure [TDl 



5. Summary 

In this article we considered the two-dimensional Ising strip of width M with surface 
fields hi and /i2 acting on the boundaries of the system. We considered only symmetric 
{hi = h2) and antisymmetric {hi = — /i2) configurations of the surface fields. 

We introduced two pseudotransition temperatures: Tj^ and Tj^. Around Tj^, 
in the antisymmetric case, the interface separating two magnetic phases moves from 
position close to one wall to the centre of the strip. At Tj^ the difference between 
the two phases disappears. The existence of these two temperatures follows from the 
properties of our solution for the free energy. We proved that and T^m have the 
same scaling properties as real transition temperatures in higher dimensions. We also 
checked scaling relations of T2 m postulated by Parry and Evans [20j. 

The major part of our analysis was concentrated on the properties of the solvation 
force. We calculated this force as a function of temperature T, surface field hi and 
strip width M. For symmetric surface fields this force is always negative (attractive). 
For strong surface fields the solvation force has a minimum above the bulk critical 
temperature of the 2D system Tc, while for small surface fields the minimum is located 
below Tc. There exists a range of surface fields for which this force has two minima. 
For antisymmetric surface fields (and hi ^ J) the solvation force changes the sign: it 
is negative for small temperatures and positive (repulsive) for high temperatures. The 
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temperature T* at which the solvation force is zero is located very close to the wetting 
temperature of the semi-infinite system T„. For large surface fields the solvation force 
has a maximum below Tc. When the surface field is decreased, this maximum crosses 
Tc. Upon further decrease of hi the maximum disappears. 

To explain these properties we proposed scaling functions in three different scaling 
regimes: at T^, at Tc, and in the case when T„ — > Tj.. 

For antisymmetric surface fields close to we found scaling in the limit M —>■ oo 
and T ^ with M {T — T^) fixed. We succeeded in finding the analytical formula 
for the scaling function and used it to explain the behaviour of the solvation force. The 
scaling function is nonuniversal, i.e. it depends on the magnitude of surface field. 

Close to Tc the scaling limit is M — > oo and T —>■ with M {T — T^) fixed. For 
both symmetric and antisymmetric surface fields the obtained scaling function is, within 
our numerical accuracy, independent of hi for hi ^ 0. We also showed analytically 
that the difference between scaling functions in both configurations of surface fields is 
independent of hi. Using properties of such obtained scaling functions we explained 
the location of maxima and minima of the solvation force around Tc for strong surface 
fields. 

The third scaling limit corresponds to hi 0, which implies ^ T^, T ^ T^. 
and M — > oo with M {T — Tc) and Mhf fixed. In this limit we calculated the scaling 
function numerically and checked that it explains the location of minima and maxima 
of the solvation force for small surface fields. 
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Appendix A. Numerical calculation of scaling function {x,y) 

In this appendix we explain methods used to calculate the scaling functions. From (!56ll 
and (|571) we get the formula 

r (x, y) = hm (x, y) , y^, (x, y) = M'fX (T (x, M) , hi (x, y, M) , M) , (A.l) 

where T (x, M) = Tc (l + x^^/M) and hi (x, y, M) = yknTjAo [x^^/mY^^. Function 
y^.j (x, y) can be calculated numerically with arbitrary numerical precision. However, 
when M is large or high precision is required, the time spend on calculation becomes 
very long. Although the limiting value 3^" (x, y) cannot be calculated exactly, it may be 
estimated in several ways. One possibility is to fix a large but finite M and assume 

y'^{x,y)^y2,ix,y). (A.2) 

This method was used in [181 with M = 200. 
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In this paper we applied the least squares method. Because the difference 
[x, y) — (x, y) depends on M and its absolute values are large for small M, 
this method cannot be used directly. 

To overcome this problem and to estimate the value of 3^" {x, y) we calculate the 
values of {x, y) for M = Mq, Mq + 1, Mq + 2, . . . Mq + m and fit the results to the 
formula 

»(-,!/)=Bo + §H-|| + ... + |^. (A.3) 

which we assume to reflect the form of leading corrections to the scaling. We take 
Bq as our estimate of 3^" {x,y)- The accuracy of this algorithm depends on values of 
parameters Mq, m and n. The larger values used the more accurate the result is; we 
used Mq = 190, m = 10 and n = 3. 

The accuracy of such obtained results may be estimated by comparing them with 
the results obtained for different values of Mq. In addition, to test our results we used 
( l58l) with A3^ calculated analytically. The obtained relative accuracy is better than 
10-1 
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